Dynamics of Bulk vs. Nanoscale WS2: Local Strain and Charging Effects 
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We measured the infrared vibrational properties of bulk and nanoparticle WS2 in order to investi- 
gate the structure-property relations in these novel materials. In addition to the symmetry-breaking 
effects of local strain, nanoparticle curvature modifies the local charging environment of the bulk 
material. Performing a charge analysis on the zy-polarized Ei« vibrational mode, we find an approx- 
imate 1.5:1 intralayer charge difference between the layered 2H material and inorganic fullerene-like 
(IF) nanoparticles. This effective charge difference may impact the solid-state lubrication properties 
of nanoscale metal dichalcogenides. 

PACS numbers: 63.22.+m, 78.20.Ci, 61.46.+W, 78.30.-j 
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INTRODUCTION 



Inorganic fullerene-like (IF) nanostructures have re- 
cently attracted attention due to their unique closed cage 
structures and outstanding solid-state lubricating behav- 
ior |H El- J us t as carbon fullerenes are nanoscale 
analogs of layered graphite, IF nanoparticles and nan- 
otubes are curved analogs of the corresponding quasi- 
two-dimensional material. Layered and nanoscale metal 
dichalcogenides are prototypes in this regard, and the 
discovery of the WS2-based family of IF nanoparticles 
(Fig.^J provides the opportunity to investigate structure- 
property relations in bulk vs. nanoscale materials. At the 
same time, the IF materials hold out the potential for 
important applications. In addition to use in recharge- 
able batteries, optical devices, and in impact-resistant 
nanocomposities, extensive mechanical properties test- 
ing demonstrates that the friction coefficient of IF-WS2 
nanoparticles is reduced up to 50% compared with the 
2H-WS2 parent compound, maintaining excellent lubri- 
cating behavior even under very high loads, ultra- high 
vacuum, and in humid conditions 
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Because of these observations, major efforts have been 
directed at understanding the connection between bulk 
and microscopic properties and exploiting the commer- 
cial promise of these novel nanomaterials. 

2H-WS2 belongs to the space group P63/mmc (-Dgfc) 
and contains two formula units per unit cell |ll|. The 
bonding is well-known to consist of strong covalent in- 
tralayer forces and weak van der Waals interactions be- 
tween adjacent layers [T2I H^ |. Each MX2 layer (M = 
group VIB metal, X = group VIA clement) contains a 
layer of metal atoms, sandwiched between two chalcogen 
layers, with each metal atom bonded to six chalcogen 
atoms in a trigonal prismatic arrangement. A group the- 
oretical analysis (Table IA~T1 Appendix) gives a total of 
18 normal modes [Tj|- The doubly degenerate E lu and 
singly degenerate A2 U vibrational modes are infrared ac- 
tive; the conjugate gerade modes are Raman active. Fig- 
ure |21 displays the relevant displacement vectors. Note 
that the infrared active xy-polarized Ei„ and z-polarized 
A2 U modes are associated with intralayer and interlayer 




FIG. 1: Transmission electron microscope images of (a) lay- 
ered 2H-WS2 and (b) a representative IF-WS2 nanoparticle 
(b). Each nanoparticle consists of a hollow core and several 
W-S-W layers. The average particle diameter is ~120 - 200 
nm for the materials of interest here. The IF nanoparticle 
shown here is slightly smaller than average. 



motion, respectively. This directional selectivity provides 
a sensitive and microscopic probe of charge and bonding 
interactions that we will employ in this work to assess 
effective total and local char ge d ifferences between the 
bulk and nanoscale materials [l5|, [l(| fl7L llSl Qjj . X-ray 
diffraction reveals that the local structure of 2H-WS2 is 
preserved within the unit cell of an IF nanoparticle with 
the exception of a 2% lattice expansion along the 2-axis 
ill. This lattice expansion is attributed to strain 
in the curved WS2 layers 0, each of which has 
a slightly different radius due to layer inhomogencitics 
(Fig. Efb)). The lattice modes of 2H-WS2 were previ- 
ously investigated using combined Raman and inelastic 
neutron scattering, demonstrating that the two-phonon 
resonance Raman effects are second-order processes in- 
volving the longitudinal acoustic mode at the K point 
of the Brillouin zone [2j| . Previous studies also indicate 
that the new Raman peaks in spectra of IF-WS2 nan- 
otubes and nanoparticles should be assigned as disorder- 
induced zone edge phonons |24[. More recent Raman 
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FIG. 2: Displacement vectors of the infrared and Raman 
active modes of 2H-WS2. The dashed line between the two 
layers represents the weak van der Waals force. 



work demonstrates that the 33 cm -1 E2 9 shearing or 
"rigid layer" mode is almost completely blocked in IF- 
WS2 due to surface strain hindering intralayer motion in 
the nanoparticles j2^. Extensive optical properties and 
STM work indicates that the indirect gap is smaller in the 
IF nanoparticles compared to that in the bulk [2(| [2?J • A 
and B exciton positions are also sensitive to confinement 
and the number of layers in the nanoparticle [24I l2fjj . 

In order to investigate structure-property relationships 
in these chemically identical but morphologically differ- 
ent metal dichalcogenides, we measured the infrared re- 
flectance spectra of both 2H- and IF-WS2 powders and 
performed a charge analysis to extract total and local ef- 
fective charge from the oscillator parameters of the major 
infrared active phonon modes. We found an approximate 
1.5:1 intralayer charge difference between the 2H- bulk 
and IF-nanoparticle materials, respectively. The trend is 
different in the interlayer direction, reflecting a slightly 
enhanced z interaction in the nanoparticles. We discuss 
how differences in both charge and strain may be con- 
nected to the macroscopic properties of these materials. 
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FIG. 3: 300 K reflectance spectra of 2H- and IF-WS2. The 
Eiu and hi u modes are infrared active in these materials, 
with xy- and z-directed polarizations, respectively. 

Elmer Lambda 900 grating spectrometer, covering the 
frequency range from 25 - 52000 cm -1 . A helium-cooled 
bolometer detector was employed in the far-infrared for 
added sensitivity. Both 0.5 and 2 cm -1 resolution were 
used in the infrared, whereas 3 nrn resolution was used in 
the optical regime. Variable temperature measurements 
were carried out with a continuous-flow helium cryostat 
and temperature controller. A Kramers-Kronig analysis 
was used to calculate the optical constants from the mea- 
sured reflectance, yielding information on the dispersive 
and lossy response of each material [3^, ■ Standard 
peak-fitting techniques were employed, where appropri- 
ate. 



III. RESULTS AND DISCUSSION 



II. EXPERIMENTAL METHODS 

IF-WS2 was prepared from its oxide precursor, WO3, 
followi ng p reviously published procedures 0, 0, |2£| |2j| 
I30L l3ll l32| . The IF nanoparticles of interest in this work 
range in size from ~120 - 200 nm in diameter. The 
particle size, shape and distribution have been studied 
by x-ray pow der diffraction |2l|, scanning tunneling mi- 
croscopy |27||. and high-resolution transmission electron 
microscopy |2j, |23, l2q |. Pressed isotropic pellets were 
prepared to investigate the dynamical properties. Bulk 
2H-WS2 was also measured for comparison (Alfa Aesar, 
99.8%). 

Near-normal infrared reflectance was measured over a 
wide frequency range using a series of spectrometers in- 
cluding a Bruker 113V Fourier transform infrared spec- 
trometer, an Equinox 55 Fourier transform instrument 
(equipped with a microscope attachment), and a Perkin 



A. Understanding Charge Localization Effects in 
2H- and IF-WS 2 

Figure displays a close-up view of the far infrared 
reflectance spectra of 2H- and IF-WS2 at 300 K. We as- 
sign the major features in the reflectance (at 356 and 437 
cm -1 ) as the Ei„ and A 2 „ modes, respectively [lH l35| . 
These modes are strikingly different in the two materi- 
als. The Ei„ mode appears damped and suppressed in 
IF-WS2 compared to the 2H- analog, whereas the A2 U 
mode is slightly more pronounced in the IF compound 
compared to that in the bulk. These differences can be 
quantified using traditional dielectric oscillator models 
and fitting techniques which, in combination with appro- 
priate models, allow us to assess the charge characteris- 
tics of the nanomaterial as compared to the bulk. 

One well-established approach for quantifying charge 
in a material involves assessment of both total and local 
effective charge HE EE E3, EE Ell As discussed by 
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TABLE I: 

Classical Oscillator Parameters and Optical Phonon 
Frequencies for 2H- and IF-WS2 







Mode 


A 2u 


Mode 


300 K 


2H-WS2 


IF-WS2 


2H-WS2 


IF-WS2 


Oscillator strength, S 


0.031 


0.014 


0.002 


0.002 


Damping constant, 7 


0.005 


0.010 


0.004 


0.006 


ujto (cm -1 ) 


356.47 


356.97 


438.32 


439.01 


ujlo (cm -1 ) 


357.10 


357.28 






ujlo-ujto (cm -1 ) 


0.63 


0.31 






£00 


9.58 


6.26 






10 K 


Oscillator strength, S 


0.040 


0.017 


0.005 


0.004 


Damping constant, 7 


0.004 


0.010 


0.006 


0.004 


ujto (cm -1 ) 


358.30 


358.61 


440.80 


440.87 


ujlo (cm -1 ) 


359.14 


359.20 






lolo-uto (cm -1 ) 


0.84 


0.59 






Coo 


8.92 


6.46 







Burstein et al. .19|. macroscopic effective charge, e^, is a 
measure of the electric moment per unit cell and contains 
contributions from both the localized charge on the ion 
sites and the charge generated throughout the unit cell. 



Therefore, ej, can be separated into a localized part, 
and a nonlocalized part, e*,, 



(1) 



Here, e* is denned as the localized moment generated 
per unit displacement of an ion. It induces a local field 
through dipole-dipole interactions which contributes to 
the reduction of the transverse optical (TO) phonon fre- 
quency. Although can be a good measure of bond 
ionicity or covalency in layered MX2 transition-metal 
dichalcogenides, e ; * gives a more appropriate representa- 
tion of bonding interactions because it quantifies charge 
on the ionic sites [l5| . 

Total macroscopic effective charge, e^, is given as 



e 



(2) 



Here, wto is the TO phonon frequency (in cm -1 ), S is 
the oscillator strength, N is the number of WS? formula 
units per unit volume, fh is the mode mass [3J3], e is 
the electronic charge, c is the speed of light (in cm/s), 
and eo is the permittivity of free space. Note that for 
comparison of two similar materials, such as 2H- and 
IF-WS2, oscillator strength and TO phonon frequency 
will be the distinguishing parameters. Localized effective 



charge ef, is determined as 



(u, 2 



LO 



■ CO, 



TO) 



me 
e 2 LN' 



(3) 



Here, oj l o is the longitudinal optical (LO) phonon fre- 
quency (in cm -1 ) and L is the Lorentz factor for a hexag- 
onal lattice I33. Precise measurement of the LO-TO 
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FIG. 4: Frequency dependence of the imaginary part of 
the dielectric function, £2(01), and the energy loss function, 
-Im(l/e(w)), of IF-WS2 at 300 K. Longitudinal optical (LO) 
and transverse optical (TO) frequencies are indicated. The 
inset displays an oscillator fit (white line) to the real (black) 
and imaginary (gray) parts of the dielectric function, £i(w), 
and € 2 (uj), for IF-WS2 at 300 K. 



splitting is key to distinguishing local charge differences 
in 2H- and IF-WS 2 . 

In order to obtain the parameters needed to extract to- 
tal effective and local charge for both 2H- and IF-WS2, 
we carried out a Kramcrs-Kronig analysis of the mea- 
sured reflectance. LO and TO phonon frequencies are 
obtained directly from the optical constants, for instance, 
as peaks in the energy loss function, — Im(l/e(u;)), and in 
the imaginary part of the frequency dependent dielectric 
function, €2(00), respectively. Figure0]shows an example 
of these optical constants and the straightforward extrac- 
tion of LO and TO phonon frequencies for IF-WS2 at 300 
K. Oscillator strength, S, is obtained by simultaneously 
fitting the real and imaginary parts of the complex di- 
electric function, e(u>) = £i{oj) + 162(00), using the three 
parameter model 
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(4) 



Here, the subscript j refers to the mode of interest, 7 
is the damping constant, and is the high frequency 
dielectric constant [1^, 13 HE 113 ■ As an example, the 
inset of Fig.0]shows an oscillator fit of the real and imag- 
inary parts of the complex dielectric function for the Ei u 
mode of IF-WS2 at 300 K. The classical oscillator param- 
eters of both 2H- and IF-WS2 are summarized in TableQ] 
These parameters, along with Eqns. [21 and [31 allow us to 
evaluate the charge characteristics. 

The small LO-TO splitting in both 2H- and IF-WS 2 
indicates that the metal-chalcogen bond is highly cova- 
lent within the layer, in agreement with previous results 
for WS2 as well as other well-known covalent compounds 
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TABLE II: 



Macroscopic and Localized Effective Charge Values for 2H- 
and IF-WS2 







Eiu 


Mode 


A 2 „ 


Mode 






2H-WS2 


IF-WS2 


2H-WS2 


IF-WS2 


300 K 


e* T /e 


0.45 


0.30 


0.15 


0.16 




e?/e 


0.20 


0.14 






10 K 


e* T /e 


0.51 


0.34 


0.22 


0.20 




e */e 


0.23 


0.20 







"e* cannot be obtained for the z-polarized interlayer A2 U 
mode due to the negative Lorentz factor |40j. 
'Error bars on the total and local charge values are ± 0.02. 



including MoS 2 , MoSe 2 , and WSe 2 d E El E E3 ■ 
In contrast, Uchida and Tanaka report large LO-TO 
splittings for several group IV transition-metal dichalco- 
genides, including lT-TiSe 2 , lT-ZrSe 2 , and lT-HfSe 2 , 
which are considered to be highly ionic materials [15J. 
Physically, the smaller LO-TO splitting of the Ei„ mode 
in IF-WS 2 (Table 1), indicates that the IF- nanoparticles 
are slightly more covalent than the parent 2H- compound. 
This result indicates that nanoparticle curvature changes 
the charge-sharing environment within the layer. 

Table[n]displays the total macroscopic and local charge 
for both 2H- and IF-WS 2 . Using the a^-polarized E ltl 
mode as a probe of charge within the layer, we find that 
4 = 0.45 for 2H-WS 2 and e* T = 0.30 for IF-WS 2 at 300 
K j3^. Thus, IF-WS 2 has approximately two thirds the 
intralayer charge as 2H-WS 2 . This approximate 1.5:1 
charge difference is replicated in the local charge num- 
bers, with e ; * decreasing proportionally in the IF com- 
pound. The intralayer charge differences summarized in 
Table ILTI can be traced to differences in the LO-TO split- 
ting and oscillator strength of the E 2 „ mode in 2H- and 
IF-WS 2 . As already mentioned, these changes are easily 
observed in the spectral data (Fig. |3J) . The results imply 
that there is in fact a significant difference in the local 
environment and chemical bonding between the 2H- and 
IF- materials and that nanoparticle curvature changes 
the charging environment within the plane. A blocked 
"rigid layer" E 2ff Raman mode is also consistent with 
these observations j25j. 

In the interlayer direction, the charge trend as char- 
acterized by the behavior of the A 2 „ mode is different 
(Table |nj. At 300 K, e* T = 0.16 for IF-WS 2 compared 
with 4 = 0.15 for 2H-WS 2 . That e* T is larger in IF-WS 2 
is indicative of slightly stronger interlayer interaction and 
enhanced charge environment in the curved nanoparticles 
compared with the bulk. As expected, the total effective 
charge within the layer (e^ from the Ei u mode) is always 
larger than that between layers (e^ from the A 2u mode), 
indicating that the majority of charge resides within the 
metal dichalcogenide layer in both materials. 



Table |n] also displays the total and local effective 
charge of 2H- and IF-WS 2 at low temperature. Although 
the exact values of total and local effective charge differ 
slightly from their 300 K values (for instance, the total 
charge extracted from analysis of the x?/-polarized E lu 
mode of 2H-WS 2 is 0.51 at 10 K and 0.45 at 300 K), 
the overall trends between the layered and nanomaterial 
remain similar to those discussed above. Macroscopic ef- 
fective charge within the layer decreases from 0.51 in the 
bulk to 0.34 in the nanomaterial (again, an approximate 
1.5:1 charge difference), and local charge makes up ap- 
proximately 50% of the total charge. The total charge 
extracted from the interlayer A 2u mode is slightly larger 
in the 2H material than in the IF nanoparticles at low 
temperature. 



B. Curvature-Induced Local Symmetry Breaking 
in IF-WSa 

Are there other manifestations of curvature in the 
nanoparticles besides the aformentioned total and local 
effective charge differences? Certainly, strain and con- 
finement have been of recent interest in both vanadium 
oxide inorganic nanotubes and silicon nanowires EJ0. 
In both cases, strain broadens the vibrational modes. 
Another effect of curvature is that the local, short range 
symmetry is formally lower than (and a subgroup of) the 
unstrained bulk. The reduction of local symmetry can 
change the selection rules, allowing formerly "infrared- 
silent" modes to become infrared active j43j,|44j. Further, 
the curvature of each metal dichalcogenide layer within 
the nanoparticle is not uniform. This inhomogeneous 
structure also results in mode dispersion. Evidence for 
these effects, while present in the 300 K spectrum of IF- 
WS 2 (Fig. is best illustrated in the low temperature 
spectral response. 

Figure displays the optical conductivity, ai(u>), of 
2H- and IF-WS 2 at 10 K 0. The low temperature 
spectral response of both materials is still dominated by 
the E ltl and A 2u modes, although because of the addi- 
tional fine structure, the nanoparticle response is clearly 
much richer and more complicated than that of the bulk 
material. We attribute the additional vibrational struc- 
ture to the formally "silent" and combination modes, ac- 
tivated (and dispersed) in the spectrum of IF-WS 2 by 
the symmetry breaking that results from the curved cage 
structure. For instance, some Raman modes are well- 
known conjugates of infrared active features [45T |. The 
E^ and conjugate pair is an example. The symme- 
try analysis and vector displacement diagrams of Verble 
and Wieting ^3 (see also Table lA~Tl Appendix) also pro- 
vide several candidates for silent mode activation, stating 
that the inactive Bi g , B 2tl , and E 2tl modes are nearly de- 
generate with several Raman- and infrared-active modes. 
Thus, the optical conductivity of IF-WS 2 likely contains 
weak contributions from all of these first-order structures, 
along with a substantial number of peaks that originate 



5 




330 350 370 390 410 430 450 

Frequency (cm 1 ) 



tween the 2H- and IF- materials. The trend is different 
in the interlayer direction (probed by the A 2ll mode) at 
300 K. Here, effective charge increases slightly in the IF 
material compared to that in the bulk. We attribute the 
stronger interlayer interaction to nanoparticle curvature. 
Elucidating the relationship between structure and effec- 
tive total and local charge in 2H- and IF-WS2 is only 
the first step in understanding the fundamental interac- 
tions underlying the phenomenal mechanical and solid 
state lubricating properties of nanoscalc metal dichalco- 
genides. Tuning MX2 bond covalency may, for instance, 
allow simultaneous exploration of macroscopic mechan- 
ical properties, the charge environment, and potential 
surface effects. 
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FIG. 5: Optical conductivity of 2H- and IF-WS2 at 10 K. 
Arrows highlight the low-temperature symmetry breaking in 
the IF material. Ref. [4^1 lists the frequencies of these fea- 
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should be noted that the exact details of this fine structure 
will likely depend on particle size and distribution within a 
sample, the number of nanoparticle layers, and the nature of 
the defects. 

from mode combinations and disorder effects. 

Group theory can also be used to investigate possible 
second order infrared active combinations. When a di- 
rect product of two fundamental modes contains A 2 „ or 
Eiu symmetry in its character, it follows that a combi- 
nation of those two modes will be infrared active. For 
example, E lu (g) E 2g = B lu + B 2u + E lu , and A 2u (g) 
Ei s = Ei u . These combinations are both infrared ac- 
tive because their reducible representation contains Ei„. 
(Taking the direct product of this state with the x, y, 
and z dipole moment operators clearly contains the to- 
tally symmetric group.) Similarly, the sum and difference 
of other fundamental vibrational frequencies can lead to 
various infrared active combinations like A 2ll (g) Ai s , Ei u 
(g) Ai g , Ei„ (g) Ei 9 , B 2tl <g> E 2g and B 2u (g) B lg . We note 
that the longitudinal acoustic mode is intimately and ex- 
tensively involved with many combination modes that 
appear in resonance Raman 23] and may be involved in 
the infrared-active combination modes of IF-WS 2 as well. 

IV. CONCLUSION 

We report the infrared vibrational properties of bulk 
and nanoparticle WS 2 in order to investigate the 
structure-property relations in these novel materials. In 
addition to the symmetry-breaking effects of local strain, 
nanoparticle curvature modifies the well-known charge 
environment of the bulk material. Using the Ei u mode 
as a local probe of charge behavior within the layer, we 
find an approximate 1.5:1 intralayer charge difference be- 
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APPENDIX: Group Theoretical Results for 2H-WS 2 

In 1970, Verble and Wieting carried out a complete 
symmetry analysis on hexagonal layered compounds with 
the goal of analyzing vibrational mode symmetries. Our 
analysis differs from that in Ref. 14] in two ways: (1) a 
C 2 correlation table accounts for the D^h 2(c) site sym- 
metry of tungsten, and (2) a ad correlation table accounts 
for the C3U 4(f) site symmetry of sulfur 0,E3- I n Ref- 
[47| refer to Table 14, p. 50, and correlation table, p. 
210. Tabic |A~T1 summarizes our results. Local strain in 
the IF-WS 2 breaks the selection rules of the 2H- "parent 
compound", as detailed in the text. 



TABLE A.l: 
Group Theoretical Analysis for 2H-WS2 



Atom 


Site Symmetry 


Irreducible Representation 


W 


D 3h 2(c) 


A2u + Bi 3 + E2 S + Eiu 


s 


G 3v 4(f) 


Alg + A2u + Big + B2 U 

+ Eig + Eiu + E2 9 + E2u 


^Total 




2A2U + 2Blg + Alg + B2u 

+ Eig + 2E2g + 2Eiu + E2U 


L Inactive 




2Bi 9 + B 2 „ + E 2 U 


r Acoustical 




A2U + Eiu 


r Raman 




2E 2 g + Alg + Eig 


^Infrared 




A2U + Eiu 
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